Quasiparticle spectra in the vicinity of a rf-wave vortex. 
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We discuss the evolution of the local quasiparticle spectral density and the related tunneling con- 
ductance measurable by the scanning tunneling microscope, as a function of distance r and angle 
9 from the vortex core in a d x 2_ y 2 superconductor. We consider the effects of electronic disorder 
and of a strongly anisotropic tunneling matrix element, and show that in real materials they will 
likely obscure the ~ 1/r asymptotic tail in the zero-bias tunneling conductance expected from the 
straightforward semiclassical analysis. We also give a prediction for the tunneling conductance 
anisotropy around the vortex core and establish a connection to the structure of the tunneling 
matrix element. 
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I. INTRODUCTION. 

While there remains almost no doubt at present that 
the hole-doped high-T c cuprate superconductors posses 
an unconventional d x 2_ y i order parameter, the micro- 
scopic origin and many phenomcnological consequences 
of this fact remain to be understood. The situation is 
most pressing in the presence of applied external mag- 
netic field where the interplay between the spatially 
varying order parameter, supercurrents, and low energy 
quasiparticles results in a great variety of novel effects. In 
the Meissner state Yip and Sauls had predicted a non- 
linear Meissner effect M manifested by an anisotropic 
component of the in-plane penetration depth linear in 
field, a unique consequence of the nodal structure of the 
d-wave order parameter. However, despite considerable 
experimental effort, this effect has not been clearly ob- 
served in cuprates j^j and, very recently, Li, Hirschfeld, 
Wolfle H argued that the effect might not be observable 
in cuprates, even in principle, due to the highly non- 
local nature of the electromagnetic response of a rf-wave 
superconductor at low temperatures. In the mixed state 
similar non-local effects have been predicted to result in 
a very rich equilibrium vortex lattice structure phase di- 
agram H), but again, no conclusive experimental confir- 
mation has yet been reported. 

Among the more successful theoretical predictions spe- 
cific to the d-wave order parameter is Volovik's predic- 
tion of a ~ T\[~H contribution to the specific heat || 
which was identified in measurements on YBa2Cu3 07 
(YBCO) single crystals [||. Although subsequent exper- 
imental investigations reported deviations from the pre- 
cise \f~H behavior |7), they were consistent with more 
general scaling relations ||||] based on the same gen- 
eral physical picture. In its simplest form the ~ Ty/~H 
behavior can be derived by observing that the super- 
current orbiting each individual vortex Doppler-shifts 
the local quasiparticle spectrum, Ek — > E^ — v s (r) • k, 
which in turn results in finite residual density of states 
N(E = 0,r) cx f s (r) cx 1/r. Integrating N(0,r) up to 
the inter- vortex distance Rh cx 1/ \H one obtains the 
total density of states per vortex and the Volovik's re- 



sult then follows on multiplying by the number of vor- 
tices n v = H/$o, with <!>o the flux quantum. We note 
that behavior consistent with the \Jli dependence of the 
residual density of states was found in numerical calcula- 
tions within the Bogoliubov-de Gennes (BdG) formalism 

Here we wish to theoretically address the surprising 
fact that the expected N(0, r) cx 1/r dependence of the 
local density of states (LDOS) has not been observed 
in the scanning tunneling spectroscopy (STS) measure- 
ments pl]-p^] , despite the fact that this state of the art 
technique definitely possesses the required spatial and 
energy resolution. Inspection of the data reveals that 
instead of a 1/r asymptotic tail at large distances the 
spectra for YBCO and Bi 2 Sr 2 CaCu 2 8 (BSCCO) re- 
cover their zero-field profiles within short distances from 
the cores on the order of several coherence lengths, be- 
yond which the spectra remain unchanged. There are 
several compelling reasons why is it desirable to resolve 
this potential conflict between the STS and thermody- 
namic measurements. The most important one has to 
do with confirming the picture of a well defined d-wave 
quasiparticle in the superconducting state of cuprates. In 
particular, since it is generally believed that the normal 
state of the cuprates is not a conventional Fermi liquid, it 
is of considerable importance to verify, in as exhaustive 
detail as possible, that the ordered state below T c , where 
the Fermi liquid scenario is believed to apply, indeed ex- 
hibits all the expected features predicted by the theory. 
When this fundamental issue has been clarified one can 
perhaps hope to tackle greater problems in the field, such 
as the nature of the non-Fermi liquid like behavior above 
T c and the origin of the pairing mechanism. 

In order to achieve our goal we first demonstrate, by 
direct comparison to the results of fully self-consistent 
BdG theory for a single vortex, that the semiclassical 
Volovik approach indeed captures the right physics of 
single-particle excitations outside the vortex core. We 
then proceed to extend this approach to the realistic case 
of quasiparticles with finite lifetimes and with strongly 
anisotropic c-axis tunneling matrix element. Our main 
result is that the above mentioned STS data away from 
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the vortex core can be understood by considering the ef- 
fect of the matrix element A4 for the electron tunneling 
along the c axis. According to the band structure con- 
siderations for tetragonal cuprates Mk exhibits the same 
anisotropy as the gap function, i.e. it vanishes linearly 
on the Fermi surface near the zone diagonals Jl4|Jl5| . In 
such a case the tunneling conductance g{V 1 r) is not sim- 
ply proportional to the temperature-broadened LDOS, 
but reflects the additional structure in Mk. Analysis of a 
case when Mk cx cos 2<p leads to a surprising conclusion 
that the the power law in the decay of the zero-bias tun- 
neling conductance <?(0, r) changes to 1/r 3 , meaning that 
it vanishes much faster than the expected 1/r decay for 
LDOS. We argue that combined with the effects of elec- 
tronic disorder, which tends to further wash out the effect 
of supercurrents on g(0, r), this mechanism is responsible 
for the absence of the 1 jr conductance tail observed in 
STS, despite the fact that LDOS itself exhibits the 1/r 
behavior. 

We further show that there exists a direct relation- 
ship between the structure of M k and the real space 
anisotropy of g(V, r) around the vortex core. In particu- 
lar we predict that the tunneling conductance along the 
circle centered at the core will exhibit fourfold anisotropy 
with maxima along the nodes and minima along the 
antinodes. In contrast, the angular dependence of LDOS 
is exactly opposite, with minima along the nodes and 
maxima along the antinodes. This prediction is directly 
testable by the STS experiment and the result should 
shed light on the structure of the tunneling matrix el- 
ement Mk which is of considerable interest in its own 
right. In the absence of disorder the ratio of minimum 
to maximum is y/2 for both g(0,r) and LDOS; disorder 
will reduce this value by adding an approximately equal 
amount of constant conductance to both the maximum 
and the minimum. 



II. TUNNELING CONDUCTANCE: GENERAL 
FORMALISM 

A. Zero field 

Tunneling conductance at bias V between a supercon- 
ductor and a normal metal is given by 

/oo 
^/'( w - e y)]T|Afk| 2 A(k, w ), (1) 
-oo k 

where / is the Fermi function and A(k, to) = 
— 2ImCJ(k, uj) is the spectral function of a superconductor 
related to the diagonal part of the full superconducting 
Greens function G(k,u>). In quasi-2D cuprates the spec- 
tral function is taken to describe electrons within a single 
ab-plane and, correspondingly, k refers to a 2D wave- 
vector. While Eq. (jl) could be written down on purely 
intuitive grounds a more detailed discussion of how one 



actually implements this dimensional reduction is given 
in the appendix. 

For the clean system the formulation based on Eq. (jl]) 
is fully equivalent to the analogous expression in terms 
BdG wave-functions but allows for a straightforward in- 
clusion of the effects of finite quasiparticle lifetime. In the 
absence of field the diagonal Greens function assumes the 
well known form |]l6| (taking ft = 1) 



0(k,w) = 



(u - iT) + e k 
( w -tT)*-e£_ A » ! 



(2) 



where £k is the normal state electron dispersion, Ak = 
Ad cos 2(f> is the d-wave gap, and T models the quasipar- 
ticle lifetime broadening that results from random disor- 
der and inelastic processes. In a d-wave superconductor, 
strictly speaking, the lifetime effects should be described 
by a frequency and wave-vector dependent self energy 
E(k, cj) whose precise structure, however, is not well un- 
derstood at present. Since we are mainly interested in 
the effect of vortices on the spectral properties at the 
lowest energies, we shall in the following ignore this com- 
plication and simply parameterize the lifetime effects by 
a constant scattering rate T. We expect this approxi- 
mation to be entirely adequate in the present context 
since the important qualitative features of the tunneling 
conductance discussed below emerge clearly in the clean 
limit r = 0. 

As alluded to in the Introduction, nontrivial depen- 
dence of the matrix element Mk on the angle <j> of the k- 
vector on the Fermi surface has measurable consequences 
for the tunneling conductance in an anisotropic supercon- 
ductor. In tetragonal cuprates band structure consider- 
ations imply strong anisotropy of the interlayer tunnel- 
ing matrix element tj_(k) oc cos2</> [ ^4[JT5| ]. One direct 
consequence of this anisotropy is the well known qual- 
itative difference between the temperature dependences 
of the in-plane and c-axis penetration depths |l4| ]. It 
is reasonable to expect that the structure of i_i_(k) will 
directly translate into similar anisotropy in the matrix el- 
ement Mk for tunneling between the superconductor and 
the STS tip. The appendix confirms this expectation by 
providing a formal derivation of A/k from the transfer 
Hamiltonian formulation of the tunneling problem. Mo- 
tivated by these considerations in the following we study 
two models: a conventional model I with A/k = M) and 
an anisotropic model II with Mk = Mi cos 2cf>, as sug- 
gested by [Q, where k = (k,cf>) in polar coordinates. In 
model I tunneling conductance is simply proportional to 
LDOS while in model II tunneling from the zone diago- 
nals is suppressed. 

It is straightforward to numerically evaluate Eq. ([l]) for 
the two models and compare g(V) to the experimental 
data. The result of the best fit to the zero-field data on 
BSCCO @ is displayed in Figure |. It is seen that, as 
pointed out previously in Ref. |17| , model II captures the 
qualitative features of the data much better than model 
I. The wide, U-shaped conductance near the zero bias 
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FIG. 1. Best fit of tunneling conductance from Eq. (|lj) to 
the zero-field experimental data on underdoped BSCCO of 
Renner et al. BM, for models I and II. The fitted parame- 
ters were Ad = 41.0 meV and V = 1.4meV for model I and 
Ad = 39.7 meV and F = 3.9meV for model II. A term linear 
in E has been added to Eq. ([[]) in order to account for the 
background conductance and the data for model I are offset 
for clarity. 



appears to be a generic feature of the c-axis tunneling 
conductance in tetragonal cuprates |T^-|2l[|, and is in- 
consistent with linearly vanishing g(V) of model I [22|. 
We therefore conclude that the available tunneling data 
are consistent with model II, as expected from the band 
structure argument presented above. 



B. Finite field - semiclassical treatment 

The effects of applied magnetic field are taken into ac- 
count by performing a semiclassical replacement 



lu - k v s (r) 



(3) 



in the Greens function of Eq. (g). Here v s (r) is the local 
superfluid velocity which in the vicinity of a single vortex 
has the form 



2mr ' 



(4) 



and is cut off exponentially at distances in excess of the 
London penetration depth A. The semiclassical approx- 
imation (||) was at the heart of the original Volovik cal- 
culation of the specific heat || and has been used exten- 
sively to compute various spectral |17 2^]2(| and trans- 
port [p7|-p9| properties of the mixed state. It appears to 
capture very well the essential physics of d-wave quasi- 
particles moving on the slowly varying background of the 
vortex lattice. 




FIG. 2. Comparison of the LDOS obtained using semiclas- 
sical approximation Eqs. (^-^) with T = (thick lines) and 
the numerical solution of the BdG equations of Ref. (thin 
lines) at indicated distances from the core. In both cases av- 
erages over the real-space angle 6 are plotted. 



We have explicitly verified that, for r > 2£, Eqs. (|l|jj) 
yield very reasonable LDOS profiles over the entire en- 
ergy range of interest when compared to the results of 
a fully self-consistent calculation within the BdG theory 
for a single <i-wave vortex |30[[5l| . Figure || illustrates the 
agreement between the two approaches. Note in partic- 
ular the excellent agreement in the low energy part of 
the spectrum. We believe that this comparison consti- 
tutes a stringent test for the validity of the semiclassical 
approximation (^) for the local quantities such as LDOS 
outside of the vortex core. As expected, however, inside 
the core the semiclassical approximation breaks down, 
as visible in the top panel of Figure |[ In the strongly 
type-II cuprates at fields well below H c i cores comprise 
only a small fraction of the total volume. Semiclassical 
approximation thus works well almost everywhere which 
explains the success of the Volovik picture in modeling 
of the mixed state. 

With the replacement (|3|) the tunneling conductance 
(|l|) becomes position dependent through the spatial de- 
pendence of the superfluid velocity. In the following we 
discuss the local tunneling conductance g(V, r) near a 
single isolated vortex in a d x 2_ y 2 superconductor. 
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III. TUNNELING SPECTRA IN THE VICINITY 
OF THE VORTEX 

A. Clean limit, zero temperature 

In the clean limit r — ► + the Doppler-shifted spectral 
function assumes a simple form 

A(k, lo)=kJ2 ( 1 - " J^) % -V + vEk) (5) 



v=±\ 



where = V e k + anc ^ V = k • v s (r). We now eval- 
uate the corresponding tunneling conductance given by 
Eq. (|l|). At T = the /' factor becomes a 5-function 
and the w-integral is trivial. The remaining sum over k 
is replaced by an integral in the usual manner. Assum- 
ing the free electron dispersion ek — k 2 /2m — ep we may 
use the (^-function in (||) to explicitly perform the inte- 
gral over the energy variable and obtain the tunneling 
conductance as a Fermi surface average of the form 



7(V,r) 

gN 



2tv 



--Re 

2?r 



2 cos 2 c 



\eV 



(6) 



Here gN is the normal state conductance, n assumes val- 
ues of and 2 for models I and II respectively and we 
have restricted k in rj to the Fermi surface. We may thus 
write 



7T £ 

i] = k F ■ v a (r) = -A d - sin( 
2 r 



(7) 



with r = (r, 8) and £ = vp/irAd the coherence length. 

For arbitrary bias V the conductance (0) must be eval- 
uated numerically. However, for \eV\ <C Ad, the integral 
is dominated by the regions close to the four nodes of 
A((j)) and can be evaluated, to an excellent approxima- 
tion, by expanding to leading order in <f> near the nodes. 
For the zero-bias conductance we thus obtain 



.9(0, r) 



9N 

where 

and 4>i = 
implied by Eq. 

g(0,r) 



2 r 



$m{6 -4>i)\, 



(8) 



(9) 



7r(2? — l)/4 are the nodal points. The integral 



|) is elementary and yields 

n+l 



sin i 



I n+l 



COS i 



lfl+1 



(10) 



where = 9 — 7r/4is the polar angle measured from a 
node, c = tt/4 ~ 0.78 and c 2 = tt 3 /16 ~ 1.94. 

According to Eq. ( p^| ) the symmetry of the tunnel- 
ing matrix element Mk has profound consequence for the 
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FIG. 3. Zero-bias tunneling conductance as a function of 
real space polar angle 9 around the vortex normalized to unity 
along the nodal direction 6 = 7r/4. 



spatial dependence of #(0, r) near the vortex core. Most 
importantly we notice that the decay with the distance 
r from the core is much more rapid in model II, where 
g(Q, r) ~ (£,/r) 3 , compared to the £/r behavior in model 
I. For instance, at r = 3£ the zero-bias conductance will 
be suppressed by a factor of 3 in model I but by a factor 
of 27 in model II. We argue that this difference is a very 
likely reason for the observed absence of 1 jr tails in STS 
measurements . 

Eq. ( |Io[ ) implies an interesting prediction for the the 
angular dependence of g(0, r) for fixed r, corresponding 
to taking a scan along a circle of the radius r centered 
at the core. This angular dependence is a result of the 
underlying k-space anisotropy of the gap function and is 
illustrated in Figure [| for models I and II. While the ra- 
tio between maximum and minimum is V2 for the both 
models, we observe that there is a qualitative difference 
between the two models in the positions of maxima and 
minima. Based on our expectation that model II pro- 
vides the correct description of tunneling into BSCCO, 
we predict that <?(0, r) should exhibit minima along the 
8 = and maxima along 7r/4 and equivalent directions. 
It would be most interesting to see if this expectation 
could be confirmed experimentally. 

It is interesting to note that numerical calculations 
within the Eilenberger formalism |32"| , |33"[ show maxima 
of zero-bias LDOS along the ir/4 directions, in contra- 
diction to the above conclusions. On the other hand, 
numerical computations within the BdG formalism |m[ , 
when averaged over a narrow range of energies near zero 
bias, show LDOS consistent with our results. 



B. Effect of temperature and finite lifetime 

It is evident that finite lifetime and temperature will 
affect the prediction (|l^) for the tunneling conductance, 
since they both lead to finite zero-bias conductance even 
in the absence of field. One expects that these effects 
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become significant at distances from the core where the 
characteristic Doppler-shift energy Ed = Ad(£/r) be- 
comes comparable to T or T. At 4.2K, which is typically 
the temperature of STS experiment, temperature broad- 
ening becomes important at r/£ > A^/4.2K which is a 
number of order 100 in cuprates. Thus, for low tem- 
perature tunneling the effects of thermal broadening are 
unimportant, except at large distances from the cores. 
Scattering rate T, on the other hand, can be a significant 
fraction of the maximum gap in cuprates, as evidenced 
by a relatively large zero-bias conductance observed ex- 
perimentally |l2|,^0|], and will therefore cause significant 
broadening at distances of several £ from the core. We 
now discuss the effect of finite lifetime T in some detail. 
For r > the spectral function becomes 



v=±l 



£ k J (oj-ri + vE^f + T 2 ' 



(11) 



Evaluation of the k-sum in Eq. (^) is now somewhat more 
involved since we no longer have a ^-function at our dis- 
posal to perform the energy integral. At T — the zero- 
bias tunneling conductance is 



1-2-K 

3 (0, r) = C / #(cos2<?!>) r 
Jo 



: / de 



(r/ + £) 2 +T 2 {^-Ef+T 2 



(12) 
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FIG. 4. Effect of finite quasiparticle lifetime V on the am- 
plitude of the zero-bias tunneling conductance as a function 
of distance r from the vortex core (for = 0). Note the 
log-log scale. The straight lines with slopes -1 and -3 mark 
the expected asymptotic behaviors for T = in model I and 
II respectively. 



where E = \J e 2 + A(<fi) 2 and C contains all the constant 
prefactors. We evaluate ([jj]) by linearizing all the func- 
tions in the integrand containing <f> around the four nodes 
of the gap function, e.g. A(</> — w 2Ad<fi, and extend- 
ing the angular integration to infinity. Under the assump- 
tion that r, \rj\ -C Ad the resulting integral can be evalu- 
ated by making use of new variables, u = \/e 2 + (2Ad(j)) 2 
and a = arctan[e/2Arf0]. The a-integration is trivial and 
we obtain 



g(0,r) 

9N 



-i 

2 = 1 



A / u X n+1 

,!,,[ a7 



(u- A dZ/ ) 2 + r 2 



( u + A dZl ) 2 + r 2 

(13) 



where z; is defined in (g) and A is a cutoff of the order 
of Ad imposed in order to assure convergence at large 
u J35j. Finally, carrying out the integral and keeping in 
mind that T, \r)\ <C A^, we obtain the leading terms in 
model I, 



g(0,r) _ 1 y 
9N 4tt ^ 



and for model II, 



2zi arctan 



7 ln(7 2 



(14) 



g(0,r) 
9n 



-t 

i=i 



2 Zl (z? 







(?) 


+ 7 







(15) 



Here we have set A = A^ and defined a dimensionlcss 
scattering rate 7 = T/A d . 

Equations ( |14| ) and (15) exhibit the correct behavior 
in the limit T — ► + when compared to the result for the 
clean system (pX|). For finite T they describe the crossover 
from the behavior close to the core (r <C rr) to 

the lifetime dominated constant zero-bias conductance 
far from the core (r ^> rr)- Inspection of Eqs. ( fTi]) and 
( |l5| ) reveals that the dependence of the crossover scale 
rr on T is more subtle than one would expect from the 
simple argument involving the Doppler-shift energy Ed 
presented above. In particular we find that 



T £/(7W) 



for model I while 



1/3 



(16) 



(17) 



for model II. In both cases lifetime effects become impor- 
tant at shorter distances from the core than one would 
expect from the naive estimate rr ~ £/7- This is illus- 
trated in Figure |] where we plot the crossover functions 
for the tunneling conductance given by Eqs. @ and (|l|) 
as a function of distance r from the core for a number 
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FIG. 5. Effect of finite quasiparticle lifetime Y on the angu- 
lar anisotropy of the zero-bias tunneling conductance. Upper 
panel: angular dependence (normalized to unity at 9 — 7r/4) 
for F/A d = 0.01 and distances r/£ = 1, 10, 20 in the or- 
der of diminishing anisotropy. Lower panel: maximum and 
minimum conductivities as a function of distance r from the 
vortex core for F/A^ = 0.01. 

of lifetimes T. We note that in model II lifetime effects 
are much more efficient in destroying the clean power law 
behavior than in model I, consistent with Eqs. (Iq) and 

The lifetime effects will also affect the anisotropy of the 
tunneling conductance around the core, since one would 
expect <?(0, r) to become isotropic at distances larger than 
rr • Figure || displays the evolution of the anisotropy as a 
function of angle 9 for increasing distance r at constant T 
(upper panel). Lower panel compares the angular max- 
imum to the minimum as a function of r. As expected, 
lifetime effects wash out the anisotropy at distances from 
the core in excess of the crossover scale rr- 

IV. SUMMARY AND CONCLUSIONS 

Our main objective was to reconcile the apparent ab- 
sence of the 1 jr tails in the local tunneling conductance 
in the vicinity of a d-wave vortex, expected on the basis 
of a simple Volovik model, with the general consensus 
that such a semiclassical picture captures the essential 
physics of the mixed state in cuprates. We argued that 
the STS measurements are likely dominated by the non- 
trivial structure of the tunneling matrix element, which 



derives from the well known band structure considera- 
tions We showed that, for the matrix element 
of the form Mk = M2 cos(2</>) (model II), in the absence 
of lifetime effects, the zero-bias tunneling conductance 
power law is modified to 1/r 3 , making it vanish much 
faster than the 1 /r tail obtained for constant Mk (model 

I)- 

We predicted a substantial angular anisotropy of the 
zero-bias tunneling conductance, with maximum to min- 
imum ratio of -\/2, and with exactly opposite arrange- 
ment of extrema in the two models. Our conclusion re- 
garding the precise form of this anisotropy, however, is 
somewhat less certain in view of the fact that it disagrees 
with the numerical results obtained within the semiclassi- 
cal Eilenberger formalism |3^j3j|] (see however Ref . |mJ ) . 
Nonetheless, our present results clearly indicate that the 
structure of the tunneling matrix element will have sig- 
nificant effect on the real-space anisotropy of g(V, r). The 
reliable information about the details of this effect can be 
obtained in a relatively straightforward manner by incor- 
porating the non-trivial matrix element into a fully self- 
consistent BdG calculation along the lines of Ref. pfj| ; 
the work on this is in progress. Within this method it 
will also be interesting to study the effect of the matrix 
element on the tunneling spectra in the vortex core, a 
problem inaccessible to the semiclassical approximation. 

Lifetime effects will cut off both the power law de- 
cay and the anisotropy at distances beyond the crossover 
length scale rr given by Eqs. (|l6|) and ([l7|), where rr is 
always shorter in model II for a given value of T. This 
means that if model II is the physically relevant one, as it 
appears to be the case for BSCCO, then the lifetime ef- 
fects may likely render the experimental detection of the 
asymptotic 1/r 3 behavior very difficult. For instance, if 
we assume T/Ad = 0.1, which is physically reasonable for 
BSCCO at low T @, then Eq. @ implies r r ~ 2.1f, 
meaning that there will be virtually no asymptotic region 
where the 1/r 3 behavior could be observed. Instead, one 
would see an onset of the constant zero-bias tunneling 
conductance just outside of the vortex core. This theoret- 
ical conclusion is in fact consistent with the STS data of 
Renner et al. [jl2|. Even in the presence of substantial life- 
time effects it should still be in principle possible to ob- 
serve the predicted angular anisotropy of g(0, r) close to 
the the core. The existing experimental data on BSCCO 
in fact show a definite hint of a fourfold anisotropy, which 
is, however, discernible only at high bias fjlj. The above 
analysis suggests that the asymptotic tails might become 
observable in the cleanest samples characterized by low 
value of the zero-bias conductance in the absence of field. 
Clearly, the effect would be difficult to discern in the ex- 
isting data on YBCO jllj which display large zero-bias 
conductance of unknown origin. 

It would thus appear that the proper inclusion of the 
anisotropy in the tunneling matrix element naturally re- 
solves the conflict between the STS and heat capacity 
measurements, since the latter is obviously insensitive to 
the structure of the matrix element. While certain de- 
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tails still await experimental verification, this seems to be 
a satisfactory tentative conclusion, especially since addi- 
tional experimental evidence emerged recently that ap- 
pears to further solidify the support for the Volovik-type 
description of the mixed state in cuprates. In particular 
measurements of the complex conductivity (using coher- 
ent terahertz spectroscopy) in the mixed state of BSCCO 
films reported by Malozzi et al. |36| were claimed to be 
consistent with the picture of the Doppler-shifted local 
quasiparticle spectra |p7ft . Although somewhat less di- 
rectly the experiments on thermal transport in the mixed 
state of cuprates |||^] also seem to support the basic 
picture of a Doppler-shifted local quasiparticle in that 
they are consistent with theoretical models based on this 
picture pEH^l . We may thus conclude that by invok- 
ing the properties of the tunneling matrix element and 
lifetime effects the existing tunneling data indeed can be 
brought to agreement with other experiments and that 
the overall picture of the mixed state in cuprates appears 
consistent with that of a well defined quasiparticle prop- 
agating on a background of slowly varying supercurrents 
in the vortex array. 
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APPENDIX: C-AXIS TUNNELING 
CONDUCTANCE IN CUPRATES 

According to the standard many-body formulation ]T^ ] 
the tunneling current at a bias V between a superconduc- 
tor and a normal metal is given by 

^A s (k,e)A N (p,e-eV) 
-co 2tt 

x[f(e)-f(e-eV)}. (Al) 



k. P 



Here As and An are spectral functions of the super- 
conductor and the normal metal respectively, T kp is the 
tunneling matrix element and e is electron charge. In the 
following we make the usual assumption that the normal 
metal in the STS tip can be described by the spectral 
function 



A N (p,e)=2ir5(e-£ p ) 



(A2) 



with a simple free electron dispersion £ p = (p^ + 
p 2 z )/2m — ep where we have resolved p into components 
parallel (jp z ) and perpendicular (p_i_) to the z axis along 
which the tunneling current flows. The 5- fun ction in A^ 
can be used to perform the e-integral in (Al) to obtain 



I(V) = 2e |T k p| 2 A 5 (k, £ P + eV)[/(£ p + eV) - /&,)]. 

k,p 

(A3) 

The expression for the tunneling matrix element can 
be derived from the transfer Hamiltonian formalism [M : 



dk z 



dp z 



£(e z )<5(k ± - PJ 



(A4) 



Here e k denotes the normal-state dispersion in the super- 
conductor, D(e z ) is the barrier transmission coefficient 
(e z is the energy of the electron inside the barrier), and 
the (5-function reflects conservation of the perpendicular 
momentum under the conditions of specular transmis- 
sion. The latter can be used to carry out t he su mmation 
over the transverse momentum p± in Eq. (A3). The re- 
maining sum over p z can be converted into an integral in 
the usual manner. The integral is then transformed by 
making a substitution u> — {p 2 z + k^)/2m — e_F + eV and 
noting that the term \d^ p /dp z \ in |Ikp| 2 is precisely the 
Jacobian of this transformation; we obtain 



/ = 2e£ 



<9e k 



D(e z ) 



2^ 



As(& U )[f(u)-f(w-eV)]. 



(A5) 



Differentiating with respect to V we finally arrive at the 
expression for the funneling conductance 



g(v) 



1 W -/V-eV)5>f k | 2 ^ s (k,c), (A6) 



' — oo 



27T' 



'■\d^/dk z \D(e z 



with |M k | 2 

Although this last expression bears similarity to Eq. 
([!]), the last term still contains summation over the 3D 
wave-vector k and it therefore reflects the full 3D band 
structure of the superconductor. In the follo wing we shall 
specialize to quasi-2D cuprates and simplify ( |A6| ) further 
by eliminating the k z summation, thereby expressing the 
tunneling conductance in terms of a6-plane properties 
only. To this end we consider the single particle dis- 
persion e k of the form deduced for tetragonal cuprates 
by Xiang and Wheatley Q : 



f k 



t z (k ± ) cos k z 



(A7) 



where e k is the a6-plane dispersion depending only on 
k i and 



t*(k±) 



(cos k T — cos k„ 



(A8) 



is the inter-plane tunneling matrix element. In the follow- 
ing we model ej^ by a free electron dispersion k^/2m-ej? 
and noting that only values of t z (k±) close to the the 
Fermi surface are important we approximate t z (k±) « 
t z cos 2 2(j>. While neither of these two approximations is 
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essential for the final result, adopting them will greatly 
simplify the ca lculations. We may now convert the k- 
sum in Eq. (A6) into an integral of the form (suppressing 
various constant prefactors), 



dk, 



k i dk 



\M k \ 2 A s (e k ,A;u>), (A9) 



where we have explicitly acknowledged the fact that 
Ag(k, (J) will only depend on the momentum variable 
through £k and A(0). We now convert the k± integral 
into an integral over the energy, noting that, according to 
our assumptions, only the component of et depends 
on k± and thus the Jacobian of this transformation is 
simply a constant. We thus arrive at the final result for 
the tunneling conductance 

g(v) = - T p-nu-eV) 

x / # / de\M{<t>)\ 2 A s {e,A;u), (A10) 



where 



\M(<t>)\ 2 = 



2e 2 m 


f°° dk z 


<9e k 


(2^) 2 , 




dk z 



D(e z ) = M| cos 2 24>, 
(All) 



and M2 is a c onstant. By converting the e and cj> integrals 
in Eq. (A10) back to a 2D k-vector sum we recover the 
expression for the tunneling conductance written down 
on purely intuitive groun ds in Eq. (Q). 

Superficially the result ( A10) appears to contradict the 
conventional wisdom that the band structure effects are 
"invisible" to tunneling because of the exact cancellation 
between the \dek/dk z \ factors in the T^ p matrix and the 
like factors resulting from the variable change k z — > e in 
the integral jh)|,[ll]]. In particular one could argue that 
had we done the £; 2 -integral in Eq. ( A9) first, such a can- 
cellation would have indeed occurred, leading to the re- 
sult with (^-independent matrix element. However, in the 
present case we note that the appropriate k z — > e variable 
change would be illegitimate because of the divergence in 
the corresponding Jacobian for k^ coinciding with one of 
the four nodes of t z ()s.±). Therefore, as pointed out pre- 
viously in connection with tunneling in cuprates ]l9fl , the 
tunneling conductance in fact is sensitive to non-trivial 
features in the band structure, such as the matrix ele- 
ment that vanishes along certain kj^ directions. 
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